Possible Indications of a Clumpy Dark Matter Halo 
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We investigate if the gamma ray halo, for which recent evidence has been found in EGRET data, 
can be explained by neutralino annihilations in a clumpy halo. We find that the measured excess 
gamma ray flux can be explained through a moderate amount of clumping in the halo. Moreover, 
the required amount of clumping implies also a measureable excess of antiprotons at low energies, 
for which there is support from recent measurements by the BESS collaboration. The predicted 
antiproton fluxes resulting from neutralino annihilations in a clumpy halo are high enough to give 
an excess over cosmic-ray produced antiprotons also at moderately high energies (above a few GeV). 
This prediction, as well as that of one or two sharp gamma lines coming from annihilations into 77 
or Z-y can be tested in upcoming space-borne experiments like AMS and GLAST. 

PACS numbers: 12.60.Jv, 14.80.Ly, 95.35.+d , 95.55.Vj, 95.85.Ry, 98.35.Gi, 98.62.Ck 



Some models for structure formation in the universe 
predict that cold dark matter clumps may have formed 
at various stages of the evolution of structure. Subse- 
quent hierarchical merging of small clumps into larger 
ones would evenually give rise to halos of galaxies. How- 
ever, if small and dense enough, some of these dark mat- 
ter clumps may have survived tidal interactions and could 
exist still today in the halos of galaxies including the 
Milky Way. A clumpy halo composed of supersymmctric 
dark matter particles (neutralinos) would reveal itself in 
conspicuous ways due to the increased annihilation rate 
of neutralinos into gamma rays, positrons, neutrinos and 
antiprotons. 

Several authors have pointed out that there is 

a discrepancy between the measured diffuse 7 ray back- 
ground in the Milky Way and the predictions of detailed 
emission models. A strong excess of photons with en- 
ergy above 500 MeV has recently been detected towards 
the galactic center || by the Energetic Gamma Ray Ex- 
periment Telescope (EGRET) on board the Compton 
Gamma Ray Observatory. Rather significant deviations 
from emission models in the same energy range seem to 
be present on a large scale as well. In a very recent anal- 
ysis, Dixon et al. jlj (hereafter Paper I) claim that in 
EGRET data there is a strong statistical evidence for a 
7 ray halo surrounding the galaxy. 

Assuming that the effect is real, some possible expla- 
nations have been addressed, such as an origin from un- 
resolved point sources (Geminga-like pulsars are feasible 
candidates), an underestimate of the inverse Compton 
emission at high latitudes and 7 rays associated with 
baryonic dark matter or WIMP annihilations. We will 
focus on the latter and investigate the intriguing possi- 



bility that the 7 ray halo might result from pair annihila- 
tions of dark matter neutralinos x, the lightest supersym- 
metric particle in the Minimal Supersymmetric Standard 
Model (MSSM) and one of the leading dark matter can- 
didates. We will analyse the compatibility of the signal 
with this hypothesis and compare with the information 
which could be extracted from other neutralino detec- 
tion methods, in particular low-energy cosmic antipro- 
tons. Recent results from the Balloon-borne Experiment 
with Solenoidal magnet Spectrometer (BESS) || indeed 
appear to show some excess of antiprotons at low energy, 
a characteristic signal of neutralino annihilations. 

We work in the framework of the Minimal Supersym- 
metric Standard Model (MSSM) as defined in Refs. |@,||. 
More details on our notation can be found in Ref. 
With some simplifying assumptions we are left with 7 
parameters which we allow to be varied within gener- 
ous bounds. The ranges for the parameters are given in 
Table Q. For each generated model, we check if it is ex- 
cluded by current accelerator constraints and if it is cos- 
mologically interesting, by which we mean models where 
0.025 < fl x h 2 < 1, i.e. where the neutralinos can make 
up most of the dark matter in our galaxy without over- 
closing the Universe. (f2 is the energy density in units of 
the critical density and the present Hubble parameter is 
lOO/i km s" 1 Mpc" 1 .) 

To get the normalization of the 7 ray flux from neu- 
tralino annihilations, we need to specify the distribution 
of dark matter in the galaxy. We first assume that the 
average dark matter density profile is spherically sym- 
metric and can be described by 
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Min 


-50000 


-50000 


1.0 





100 


-3 -3 


Max 


50000 


50000 


60.0 


10000 


30000 


3 3 



TABLE I. The ranges of parameter values used in our scans 
of the MSSM parameter space. Note that several special scans 
aimed at interesting regions of the parameter space has been 
performed. In total we have generated about 55000 models 
that are not excluded by accelerator searches. 



where po is the local halo density, Rq our galactocentric 
distance and a some length scale. The functional form 
of p for the Milky Way cannot be fully determined from 
observational data 10 and several choices of the parame- 
ters are acceptable. A hint on this choice may come from 
N-body simulations of hicrachical clustering in cold dark 
matter cosmologies which favour profiles that are singular 
towards the galactic centre. We will consider the Navarro 
et al. profile jll| which is cuspy and has 7 = 1 (it scales 



as 1 jr at large distances) and, for comparison, the mod- 
ified isothermal distribution, (a,(3,j) = (2,2,0), exten- 
sively used in dark matter detection computations. We 
have checked that the Kravtsov et al. profile jl^] which 
is mildly singular with 7 ~ 0.2, also gives acceptable fits. 

To make the analysis compatible with the existence 
of clumps in the halo, we consider the possibility that 
a fraction / of the total dark matter, rather than be- 
ing smoothly distributed in the halo, is concentrated in 
clumps. Simulations of structure formation in the early 
Universe do not yet have the dynamical range to give 
predictions for the size and density distribution of small 
mass clumps (we focus here on clumps of less than around 
10 6 solar masses which avoid the problem of unaccept- 
ably heating the disk jul). The formation of clumps 
on all scales is however a generic feature of cold dark 
matter models which have power on all length scales. If 
self-similarity is a guide, galaxy halos may form hierar- 
chically in a similar way to that of cluster halos (see e.g. 
|0 ) . The main effects of a clumpy halo can be sketched 
by simply assuming that clumps of given mass have an 
enhanced and roughly constant dark matter density p c \. 
We introduce the parameter S as the contrast between 
the dark matter density in clumps and the local halo 
density which is about 0.3 GeV cm~ 3 . We further as- 
sume that the average clump mass distribution follows 
the same distribution as the smooth component, Eq. (^). 

Pair annihilations of neutralinos in the dark matter 
halo can produce photons which are either monochro- 
matic or with a continuum energy spectrum. The 
monochromatic 7s arise from the loop-induced S-wave 
annihilation into the 77 or Zj final states; the phe- 
nomenology of these processes has been studied recently 
in Ref. |l5| . The continuum contribution is on the other 
hand mainly obtained from pions produced in jets; these 
photons are expected to be much more numerous but 



lower in energy than the monochromatic ones |T(|. To 
model the fragmentation process and extract informa- 
tion on the number and energy spectrum of the photons 
produced we have used the Lund Monte Carlo Pythia 
6.115@. 

Consider a detector with an angular acceptance Af2 
pointing in a direction which forms an angle ip with re- 
spect to the galactic centre. The integrated 7 ray flux 
above an energy threshold E t h is given by 



$ 7 (£ th , Ml, tP) ~ 1.87-10- 8 <S(£ t h) 
( J(if>)) (AO) cuT 2 s" 1 sr" 1 . 
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In this formula we have defined a particle physics depen- 
dent term 
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where M x is the neutralino mass, F are the allowed final 
states and for each of these vgf is the annihilation rate, 
while dN^/dE is the differential distribution of produced 
photons. The dependence of the flux on the dark matter 
distribution is contained in the factor ( J (ip) ) (Afi). In 
a scenario with many unresolved clumps, it is possible 
to show that the contribution from the clumps is given 

by pl 
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in contrast to the smooth case jl5| where the dependence 
is quadratic in the density p. The relative strength of the 
smooth and the clumped components is determined both 
by the halo profile and by the parameter fS, the product 
of the halo fraction in clumps and their overdensity. 

In Fig. [I] we show the values of 5(1 GeV) versus the 
neutralino mass for our set of supersymmetric models. 
Our results show a very large dispersion for the possible 
values of 5, about seven orders of magnitude; the highest 
S are for low neutralino masses, M x ~ 40 — 60 GeV, while 
for heavier neutralinos there are both the 1 /M ^ suppres- 
sion and lower photon production rates. The highest val- 
ues of the 7 ray flux are given by gaugino-like neutralinos 
and mixed neutralinos, whereas S is generally much lower 
for higgsino-like neutralinos, at least in the low mass 
range (the opposite trend was noticed for the monochro- 
matic photon flux, see Ref. H| ) . For the highest possible 
values of the flux S scales as l/(fl x h 2 ) which reflects the 
fact that to a first approximation fl x h 2 oc l/(va). There- 
fore, the maximal S depends crucially on the minimal 
Q, x h 2 (which we take to be 0.025) that we judge to be 
cosmologically interesting. Current estimates of Q, and h 
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a) Navarro et al.: 
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FIG. 1. The value of S versus the neutralino mass. 

indicate (with large errors) flh 2 ~ 0.1, of which baryons 
may contribute around 0.02 — 0.03 p9| . 

We focus now on the result reported in Paper I. The 
value of the residual flux at high latitudes shown in Fig. 3, 
Paper I is about 10 -6 ph cm~ 2 s _1 sr" 1 . Even pick- 
ing the MSSM model in our set which gives the highest 
<S(1 GeV), it is not possible to reproduce the qualitative 
features in that Figure in a smooth halo scenario with 
any of the three halo models considered and any of the 
allowed values for length scale a and local halo density 
po ■ On the other hand a many unresolved clump scenario 
may be compatible with the results of Paper I. We will 
show this through two examples. 

Let us first consider an example in the low mass range, 
a MSSM model which has M x = 76 GeV, n x h 2 = 0.03, 
Z g = 0.18 and which gives 5(1 GeV) = 1.025. We show 
in Fig. H for this model the values of ( J (ip) } needed to 
fit the results in Paper I. We are considering a roughly 
spherical 7 ray halo and deriving the angular distribution 
from Fig. 3, Paper I at zero longitude. As can be seen, 
the present EGRET data can be quite well reproduced 
in our approach. Although we are not in the position 
of discriminating among different halo profiles, one may 
notice that in the case of the Navarro et al. profile, a 
sharp enhancement of the 7 ray flux is predicted towards 
the galactic centre, due to the singularity of the smooth 
density profile. This is a feature which may be searched 
for in the EGRET data (some indications of results going 
in that direction were given in Ref. ||). 

We find that the dark matter signal, which causes an 
excess of photons mainly in the energy range between a 
few GeV and M x , may explain the discrepancy between 



FIG. 2. The value of J(ip) for a) the Navarro et al. profile 
and b) the isothermal sphere. The clumpyness and halo model 
parameters to fit the EGRET data, shown as filled circles, are 
given in the figure. 

the 2.75 power law fall off which is expected according to 
models of diffuse galactic background and the behaviour 
that has been mapped by EGRET up to 20 GeV. Future 
experiments [ p0|j2l] ] will be able to measure the diffuse 
background at much higher energies and eventually de- 
tect if there is a break in the energy spectrum at about 
the neutralino mass. 

For the MSSM model we have chosen, the gamma lines 
from annihilation into 77 and Z'y are well above the 
background, and their detection might be possible, es- 
pecially if there is an enhancement of the dark matter 
density towards the galactic centre. The detection of 
lines, which have no plausible astrophysical background, 
seems to be the natural way to show conclusively whether 
the 7 ray halo originates from dark matter neutralino an- 
nihilations. On the other hand, the correlation between 
continuum and line signal strengths is very weak and the 
absence of the lines would by no mean imply that there 
cannot be a continuum signal. 

Given the value of fS needed, we have to worry about 
other dark matter searches and make sure that these 
models are not excluded already. It turns out that the 
signals are below present detection limits, except for the 
antiproton signal which we now discuss. 

We have computed the antiproton signal in much the 
same way as the continuous gammas by using the Lund 
Monte Carlo Pythia 6.115. We have applied for the 
propagation the leaky box approximation with the energy 
dependent escape time given in Ref. p2| and used the 
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solar modulation model of Ref . |2^] . 

For the example considered above the antiproton flux 
at 0.4 GeV in a clumpy halo with fS = 18 is cj)p — 1.6 x 
10~ 5 cm~ 2 s _1 sr^ 1 GeV -1 which should be compared 
with the flux measured by BESS § 4> p = 1A±% x 1CT 6 
cm -2 s _1 sr _1 GeV -1 at the same energy. It may be 
tempting to conclude that this model is already excluded 
since it gives a too high antiproton flux, but one has to 
keep in mind the big uncertainties involved, mainly in 
the antiproton propagation. In particular, it is not clear 
how large a fraction of antiprotons generated in the halo 
(i.e. outside the galactic disk) can penetrate the wind of 
cosmic rays leaving the disk fzj}| . It is interesting that 
the antiproton flux is within an order of magnitude of 
the reported BESS flux, which shows some of the charac- 
teristics (pile-up at low energy) expected for neutralino- 
induced antiprotons ||. 

Since we have found in our calculations that the an- 
tiproton flux strongly correlates with the continuous 
gamma flux for a given supersymmetric model, it seems 
impossible to reduce the antiproton flux maintaining high 
continuous gamma flux at sub-100 GeV neutralino mass. 
If the overproduction of antiprotons seems uncomfort- 
ably high, it is however possible to resolve this by going 
to higher neutralino masses. 

We thus choose as our second example an MSSM model 
which has a large mass, M x = 503 GeV, fl x h 2 = 0.03, 
Z g = 0.04 and which gives 5(1 GeV) = 0.05 ph. The 
necessary rescaling for this model is fS = 427 for which 
the antiproton flux at 0.4 GeV is 4>p = 1.7 x 10~ 6 cm~ 2 
s _1 sr _1 GeV -1 , i.e. within the la error bars of the BESS 
measurement. For these higher mass models, there may 
be a problem in explaining a high gamma ray flux in 
the lower energy interval 0.3-1.0 GeV (this interval is 
however subject to larger uncertainties). An intermediate 
mass model in which the antiproton flux limit by BESS 
is mildly violated might of course be considered. More 
data is obviously needed to make firmer statements. 

We remark that for this particular model, the spin- 
independent cross section on nucleons is 0.15 • 10~ 4 pb, 
very close to the limit given by the most sensitive Nal 
experiment |25[ . For other models giving similar gamma 
and antiproton rates, however, the direct detection rates 
are much lower. 

It is interesting to note that for the high-mass neu- 
tralino, also a high-energy excess of antiprotons is poten- 
tially measurable. This is illustrated in Fig. |§|, where a 
compilation of present data ( B and references therein) 
is shown together with the predictions of cosmic-ray in- 
duced background (the mid-range of the predictions of 
p6[), and the flux in our second example. As can be 
seen, present data are not yet conclusive. However, an 
interesting feature of the high-mass neutralino result is 
that the maximum of the antiproton flux is shifted to- 
wards higher energies by 1-2 GeV compared with the low 
mass case. Also, the fall-off with energy above the peak 
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FIG. 3. The energy spectrum of antiprotons for the second 
model example (see text for details). 



is considerably slower than for the background. Indeed, 
our second model fits quite nicely this higher-energy part 
of the present data. It should be noted that at these 
energies, the effects of galactic and solar wind modula- 
tion is less severe than at sub-GeV energies, making the 
predictions more trustworthy. These features should def- 
initely be investigated in the upcoming antiproton mea- 
surements p7|p0[ . 

Although the interpretation of the measured excess in 
cosmic gamma rays and antiprotons in terms of neu- 
tralino annihilation contains elements of speculation at 
the present time, it is reassuring that upcoming experi- 
ments will be in a position to more firmly confirm or rule 
out this hypothesis. For instance, the proposed Gamma- 
ray Large Area Space Telescope (GLAST) |2l) will have 
spectral and angular resolution enough to search for 
gamma ray lines in the direction of the center of the 
galaxy. Also, if the explanation lies in a clumpy halo, 
a large-exposure experiment like GLAST eventually may 
resolve individual large clumps as bright gamma-ray 
spots on the sky. The antiproton spectrum will soon be 
measured with higher accuracy as well in the Alpha Mag- 
netic Spectrometer (AMS) @ and PAMELA |§7| ex- 
periments. Also, upcoming direct detection experiments 
may in favourable cases be sensitive to these dark mat- 
ter candidates. Finally, improved iV-body simulations 
of structure formation in cold dark matter models may 
give a better assessment of the credibility of clumpy halo 
models. 
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